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We highlight a viable mechanism leading to the formation of dark energy structures on sub-horizon 
cosmological scales, starting from linear perturbations in scalar-tensor cosmologies. We show that the 
coupling of the dark energy scalar field, or Quintessence, to the Ricci scalar induces a "dragging" 
of its density perturbations through the general relativistic gravitational potentials. We discuss, in 

■ particular, how this process forces dark energy to behave as a pressureless component if the cosmic 

' evolution is dominated by non-relativistic matter. 

, This property is also analyzed in terms of the effective sound speed of the dark energy, which 

^Nj ' correspondingly approaches the behavior of the dominant cosmological component, being effectively 

» I I vanishing after matter-radiation equality. 

CL(. To illustrate this effect, we consider Extended Quintessence scenarios involving a quadratic coupling 

' between the field and the Ricci scalar. We show that Quintessence density perturbations reach non- 

linearity at scales and redshifts relevant for the structure formation process, respecting all the existing 
constraints on scalar-tensor theories of Gravity. 

This study opens new perspectives on the standard picture of structure formation in dark energy 
cosmologies, since the Quintessence field itself, if non-minimally coupled to Gravity, may undergo clus- 
1^-^ ' tering processes, eventually forming density perturbations exiting from the linear regime. A non-linear 

' approach is then required to further investigate the evolution of these structures, and in particular 

, their role in the dark haloes surrounding galaxies and clusters. 
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CN . I. INTRODUCTION 

o _ 

The role of vacuum energy in cosmology is receiving a renewed interest since at least three cosmological observables 
CIh' indicated that a relevant fraction of the energy density in the Universe is presently in the form of a sort of vacuum 
' ' energy, which is commonly known as dark energy. Type la Supernovae observations suggest that cosmic expansion is 
2 , accelerated p],^ ; recently it hag been also noticed that data indicate that acceleration is a relatively recent occurrence 
> • in the cosmological evolution y]. Moreover, present data on Cosmic Microwave Background (CMB) anisotropics favor 
^ ' a total energy density which is very close to the critical value and which is made at roughly 70% by a vacuum energy 
component [0 . Finally, large scale structure observations suggest a universe with a sub-critical matter density |^ . 

These evidences represent a great stimulus for theoretical work, because understanding why the vacuum energy is 
at the level of the critical density today or less represents one of the main mysteries in modern fundamental physics 
At a quantum level, a "fine-tuning" mechanism is needed to explain why vacuum energy is so low with respect to 
?H any natural scale of vacuum expectations of fundamental fields: this discrepancy reaches 120 orders of magnitude if 
5y , the present cosmological critical density is compared with the Planck scale. Moreover, if the mentioned cosmological 
observables are interpreted correctly, it is necessary to justify the "coincidence" with which vacuum energy is starting 
to dominate the cosmic expansion right now. IDescribing dark energy as a scalar field (f), or Quintessence, first 
considered in 7M. has the attractive feature to alleviate these problems, at least at classical level; "tracking" 
and "scaling" ? OT] solutions existing for Quintessence with inverse power law and exponential potentials respectively, 
allow dark energy density to be at the level of other cosmological components in the very early Universe. /c-Essence 
models involving a generalized form of the kinetic energy of the scalar field can provide a mechanism to justify the 
present level of dark energy [Tl|JT^ . 

Dark energy cosmologies have oeen considered in the general framework of scalar-tensor theories of Gravity. These 
"Extended Quintessence" scenarios (hereafter EQ, while we shall refer to minimally coupled Quintessence as Q 
models), in which the scalar field responsible for cosmic acceleration possesses an explicit coupling with the Ricci 
scalar, have been studied under several perspectives [14-09], including a detailed study of tracking trajectories and of 
their effects on cosmological perturbations pO| ; in particular, both for Q and for EQ models, characteristic si gna tures 
have been accurately predicted on the CMB spectrum of anisotropics, and compared with existing data, see ||21| and 
references therein. 

Elevating the cosmological vacuum energy to a dynamical role through its representation as a scalar field introduces 
the natural question about its relation with the other main dark component which is currently under study in 
cosmology, i.e. the Cold Dark Matter (CDM). The latter is currently thought to be made of non-relativistic particles, 
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possibly generated during the process of supersymmetry breaking in the early Universe, see e.g. p^ , which are 
supposed to form the halos around galaxies via gravitational collapse of primordial linear perturbations. It is natural 
to ask if the formation of matter clumps in the Universe can have some effect on dark energy, or even if dark energy 
clumps can form starting from primordial linear perturbations. In the past this intriguing issue has been faced under 
different points of view, either by investigating the structure and stability of non-linear spherically symmetric scalar 
field overdensitiesCq] , or by describing dark matter and dark energy as two different smoothly distributed classical 
fields [^^; in j2M7clumps of an extremely light classical scalar field have been proposed as a candidate for Dark 
matter m galacticnalos. In |2^, the formation of matter structures has been faced in a background filled by matter 
and a homogeneous Quintessence, which still doesn't undergo structure formation processes. 

In any case, the occurrence of eventual da,rk energy clumps arising from perturbations in linear regime did not receive 
any explanation. Even more, in Ref. P8[ , an effective fluid with negative equation of state has been considered as 
"generalized dark matter" in the general context of cosmological linear perturbation theory, and the main consequences 
for structure formation theory have been obtained. In particular, it has been shown that perturbations of a minimally 
coupled scalar field playing the role of Quintessence behave as scalar radiation on sub-horizon cosmological scales, 
relativistically dissipating scalar field density fluctuations, so that ordinary Quintessence rapidly becomes a smooth 
component. In the present work we investigate the sub- horizon dark energy perturbation behavior in scalar-tensor 
cosmologies, and we show that the conclusions can be much different. In particular, our aim is to give an answer to 
the following question: in which conditions it is possible to form growing dark energy density perturbations on sub- 
horizon scales, in a reference-frame independent manner? We consider this problem in the general context of linear 
perturbation theory in scalar-tensor cosmology. We study the perturbation properties of the non-minimally coupled 
dark energy field, focusing on the influence that metric fluctuations can have on its sub-horizon behavior; we estabilish 
in particular in which conditions such influence is effectively dominant, resulting in a "gravitational dragging" of the 
dark energy itself. Moreover, we give a worked example of this phenomenology considering a typiecal Extended 
Quintessence model and giving numerical results on the sub-horizon behavior of its energy density fluctuations at 
redshifts relevant for structure formation. 

The results presented here are complementary with respect to our previous works on the same topic [?] : indeed, in 
p3[ we wrote the basic perturbations equations, analysing the main cosmological features of EQ models. In [ pO| , we 
focused on the impact of tracking EQ trajectories on the CMB anisotropics and on the matter power spectrurn. Here 
we deal with the clustering properties of the Dark energy itself, therefore completing the picture of linear perturbation 
theory in EQ cosmologies; this will require working in a setting in which the stress-energy tensor of the Quintessence 
field is conserved, as opposite to the formalism adopted in the previous works. 

The paper is organized as follows. In Section O we review the general formalism for scalar-tensor cosmologies. In 
Section [II we study the motion equations, both lor background and perturbations, giving emphasis to the role of the 
non-minimal coupling in the dynamics of the scalar field density fluctuations. The resulting dark energy clustering 
properties are illustrated in Section where we expose some examples of these effects as the result of numerical 
integrations in typical EQ cosmologies. In Section M w e g ive an equivalent interpretation of these results in terms of 
the effective dark energy sound speed. Finally, Section VI contains the concluding remarks. 



II. SCALAR-TENSOR COSMOLOGIES 



In this Section we give general definitions and formalism for describing general scalar-tensor cosmologies, both for 
background and linear perturbations. We follow, as much as possible, the notation adopted in original works p9| , ^ . 
Scalar-tensor theories of gravity are generally represented by the action 
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where R is the Ricci scalar and (/) is a scalar field which is supposed to be coupled only with Gravity through the 
function R), while the functions w(0) and V{(j)) specify the kinetic and potential scalar field energies, respectively; 
the Lagrangian L fluid includes all the components but 0, and the constant k plays the role of the "bare" gravitational 
constant G* , which in scalar-tensor theories can differ from the Newton's constant G as it is measured by Cavendish- 
type experiments jsi]] ; without loss of generality, we choose the relation between k and G, defined in jsi]] to be 
K = SttG*. We also pose the light velocity c equal to 1. Einstein equations from the general action (|^) can be written 
in the familiar form 
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with the stress-energy tensor T^J being made of the scalar field and the other components, indicated with the 
subscript fluid: 
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As a consequence of the contracted Bianchi identities T^"*"' is conserved; moreover, T/^"*'' and T^4(l>] are separately 
conserved since no explicit coupling is assumed between fluid and 0: 

V%^i""^ = V''T^,[(/.] = 0. (4) 

By defining 



the conserved scalar field contribution assumes the form 



Vg^, + ^^^-j^g^u + V^V.F - g^.OF +{^-F] G^, , (6) 



where we can recognize the origin of the different terms composing the scalar field stress-energy tensor (^: the 
minimal-coupling (including a "kinetic" and a "potential" part), the non-minimal coupling including /, F and i?, and 
the gravitational term, proportional to (k^^ — -F), contaning the Einstein tensor itself. We can define them as 

T;^,= [0] = u[V^c^V,c^ - ^3m.VVVa0] - Vg^, , (7) 
T^ri't^] = ^-^^^^9,u + V^V.F - g,,UF , (8) 

T3,7^m^{^^-F^G,. . (9) 

It is relevant to note that, as extensively discussed in the gravitational term may also be taken to the left hand 
side of eq. (^) : 

FG^j,^ = f^°*°' = T^''"^ + T'""[0] + T""^[(/.] (10) 

With such approach, one would be left with a not conserved total stress-energy tensor T^°*°', which would differ 

from (H) because the absence of the merely "gravitational" term in (^. Including the gravitational sector in the 
stress-energy tensor of the scalar field is not only a dichotomy, as we will discuss in this paper. First of all, in 
typical non- minimal coupling models in which _F is plus a function depending explicitely on 0, the latter function 
describes the energy transfer between field and metric, acting in particular at a quantum level, see e.g. [^3| . Second, 
the Bianchi identities allow us to write down conserved quantities for the scalar field. Third, the linear perturbation 
theory describes the evolution of small perturbation in the scalar field energy density: if the latter is drawn from 
a conserved stress-energy tensor, perturbations in the scalar field energy density exhibit a behavior which is hidden 
with the approach (^. 

Assuming a Friedmann Robertson Walker (FRW) background, the metric tensor assumes the form 

5^1/ = a^{-n){li^u + h^,^) , (11) 

where rj is the conformal time, 0(77) is the scale factor with conformal time derivative a expressed through the Hubble 
parameter Ti, — a/ a and 7^^ = diag[—l^ (1 — ifr^)"^, r^, sin^ 9] is the background metric with spatial curvature K 
in spherical coordinates {ri,r,6,(l)); h^^, <^ 1 represents linear metric cosmological perturbations, and is conveniently 
expressed in the Fourier space: for scalar perturbations, indicating with Y the solution of the Laplace equation 
ysjysy _ — j/cpy, whcrc mcaus covariant derivative with respect to the spatial metric jij, the amplitude at wave 
vector k of the most general scalar metric perturbation, see e.g. [ p9| , can be written as 

hoo = ~2AY , hoj - -BYj , /ly - 2HlYj,, + 2HtY,, , (12) 

where A,B,Hl and Ht represent the amplitude in the Fourier space at wave vector fc; Yj and the traceless Yij, 
with i,j = 1,2,3, are defined as Yj = -(l/A:)V|y, = {l/P)VfV''Y + j.^Y/S. Note that, as customary, we 
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intentionally do not write the argument k explicitly in the amplitude of the perturbation quantities in the Fourier 
space. As it is well known [ p9[ , a gauge freedom exists because of the linearization of the problem, so that two of 
the four quantities in (n3) can be set to zero, or, equivalently, two independent gauge invariant combinations can be 
built out of A,B,Hl ana Ht- Correspondingly, the stress energy tensor T^^ relative to any cosmological component 
X splits in a background component Tj^^ = diag[— p^, Px , Px , Px] and perturbations (JT^^, represented as 

STl^x = -Px 5x Y , ST"^, = {px + Px){vx - B)Yj , 5T^^ = Px{ttlx5) + tttxY;) , (13) 

where in particular 5x represent the density contrast fluctuation at wave vector k. It is also useful to introduce the 

fluctuations in the expectation value of the scalar field at wave vector fc, which will be indicated as d4)Y . 

Einstein and conservation equations (pP) split into two separate sets describing the evolution of background and 
perturbations. In the next Section we'll write them explicitly, focusing on the role of the different quantities composing 
the scalar field stress energy tensor T^i, [0] defined in (||) . 



III. GRAVITATIONAL DRAGGING 

Let us then concentrate on the conserved tensor (^, first considering background quantities and their evolution 
equations. Conservation laws (jj) for the unperturbed scalar field reduce to 

~2>H{l + w^)p^ , (14) 

having defined the dark energy field equation of state as = p^/ p^. The background evolution equations will be 
completely determined by the Friedmann-Robertson- Walker equations 

9 2 

'H^ ^ — {p fluid + Pcp) + K , H ^ ~~ ^Pn^id. + P0 + "ip fluid + 3p0) , (15) 

and by the conservation equations for each component x: Px = —3Ti,{px As it can be easily seen from (|^), the 

expressions for the scalar field energy density and pressure which satisfy dl4) are given by 



^^="£^-^(^) — 2d^+^{^^+n — — u-^j • (1^) 

It is useful to mention that these expressions combine in the continuity equation (|lj) to give the Klein-Gordon 
equation: 

4> + 2n^+^ - + 2aV^^ = . (18) 

Before considering perturbed quantities, it is relevant to comment briefly on the role of the gravitational component 
of the scalar field stress-energy tensor (Q), since as we'll see in a moment the same arguments hold for perturbations. 
As first noted in | |l8| , p^ , under conditions in which F differs from , even by a small amount due to a non-zero value 
of (p, the gravitationaTterm appearing in the expression (fiq) switches on, feeding the scalar field energy density itself 
with a term proportional to the square of the real time Hubble parameter H — Ti/a, which in turn is proportional to 
the total cosmological energy density through the Einstein equations. Since the latter is made of matter and radiation 
scaling as and 1/a^ respectively, it is straightforward that at sufficiently early times the gravitational term 

dominates the dynamics of p^. As we'll see, this process, which can be meaningfully named "gravitational dragging" , 
is also very important for the dynamics of the dark energy perturbations. 

Equations V^(5r^[^] = and V^(5Tj'[(/)] = correspond respectively to the continuity and Euler equations 

kv^ + 3Hl + m-^^T^ = , (19) 
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{v^ ^B)+H{v^- B) (1 - icl^) - k— + 

where we have defined the scalar field entropy perturbation 
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TTTtp - kA ^ , 



(20) 



(21) 



and its sound velocity 



(22) 



Equations (|I9£Q) hold formally for any cosmological component. As for the background, they combine in the per- 
turbed Klein-Gordon equation 



(2n + ^(l)]6(f> 
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A + <i) (SHA - 3Hl - ks) + -^J^c^rSR 
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(23) 



with the variation of the Ricci scalar given by 
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6R ■ 



-2 (sua - 3Hl - kB^ Y - -^n (^HA - 3Hl ~ kB^ Y- 



(24) 



In order to gain insight into the behavior of the scalar field perturbations specifically, let us write explicitly the 
formal solutions to the above equations. The variation of the stress-energy tensor (|^) yields contributions which we 
classify in mc, nmc and grav as we did in equations (0,||j9|). Let us start from the energy density perturbations: 
5ro°[(/)] = -p^S^Y = 6Tl^""'[(t)] + (5ro""™=[0] -I- JTq ^'■""[0]; the different contributions are given by 
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m^A - HkB - SHHl - {k^ ~3K) [Hl + ^ 



Y 



(27) 



where we have defined Kg = —A + H ^B/S + H ^Hl; note also that in general 5F — F^^64> + F^rSR. The momentum 
perturbation STj[(j)] is composed by 



(28) 
(29) 
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Let us consider now 5T^. In general, it will have both trace and traceless components, ttl^ and nxtf, respectively, as 
in equation ([T^). The different contributions are given by 



Y51 



(31) 
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2 .. ■ k^-3K / Hi 
-kUB ~Hl- 2nHL [Hl + -j- 



YS] + 



-k^A -k{B + HB]+HT + n{2HT - kB) - k^ [ Hl + — 



Y! 



(33) 



where tt^^ and ttj^^ are given by the terms proportional to (5* and Y^', respectively. It is worth to note that an 
interesting feature of non- minimally coupled scalar fields is the presence of the gauge invariant anisotropic stress ttti: 
as shown in |p8| and [|34| , stress perturbations have a role in the structure formation; we will return to this in Sec.|y|. 

Although the expressions above appear complicated, it is quite simple to highlight the point we are interested 
in. In the quantities ( p7| , ^ , p3| ) , the terms multiplying (1/k — F) are JGq, (JG^, (JG* , respectively. Focusing on the 
gravitational part of the scalar field energy density perturbation, as we noted above for the case of the background 
quantities, if F differs from 1/k the scalar field density perturbation is fed by the total density fluctuation, since 
6Gq = k6Tq. Therefore, if this term dominates over the others {mCjUmc), the gravitational dragging is active on 
the density perturbations and forces the scalar field density fluctuations to behave as the dominant cosmological 
component. 

As we'll see in the next Section, this process becomes crucially important in dark energy cosmologies, where the 
scalar field plays an important role in the cosmic evolution, determining the cosmic acceleration today. In the next 
Section we will give a worked example of this issue. We integrate Einstein and conservation equations to get the 
cosmological evolution in a typical Extended Quintessence scenario where 1/k — F (x (f)^ , focusing on the behavior of 
the dark energy density fluctuations at the relevant redshifts for structure formation. 
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FIG. 1. Redshift scaling of cosmological components in tracking EQ scenario (see text): radiation (dashed), matter (dotted), 
total dark energy as from eg. (p^ (solid), kinetic and potential dark energy contributions (first two terms in eq.(p^) (dotted 
dashed). 



IV. DARK ENERGY CLUSTERING 



Let us focus on the non-minimally coupled scalar field fluctuations; combining (E5|, 26|,^7|) and the expression for the 
background energy density (^6|), one gets the following expression for the scalar fiSdenergy density fluctuation: 



where 



0^ = 



+ 2aW + RF- J/k - mP - 2a?{l/K - F)Gg ' 



(34) 



(35) 



srnrac 
0^ = 



QAUF - enSF ~ a^{f^4,lK)5(t> - {-a^R + 2P)SF - QUJCgF 



+ 2a'^V + RF- f/K - enP ~ 2a2(l/K - F)G\ 



(36) 



-6(5F7^2_2a2(i/K_F)(5Gg 



+ 2a^V + RF- f/n - mF - 2a'^{l/K - F)Gl 



(37) 



Again, we point out that this is precisely the form obtained perturbing the field energy density, whenever the latter 
is drawn from a conserved energy-momentum tensor: only in this case, we are allowed to use eq. (|l9|) for the field 
energy density evolution. Most importantly, the use of conserved quantities allows a more direct interpretation of 
the interchange between different species. The gauge invariant total density perturbation A and the gravitational 
potential defined by 



Px5x 
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Px)v^ 



- {P + P)B 



(38) 
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are related through the Einstein equation 



3 A: \ k 



^ = WP^- (40) 



In such a way, since A sums up perturbations in aU the fluid components, a "potential well" may be generated by 
each of them (in particular, by a perturbation in the scalar field energy density), affecting the behavior of density 
perturbations in all the other species (in particular, matter perturbations). Viceversa, perturbations in the matter 
component will perturb the gravitational potential to drive the scalar field energy density perturbations: such a kind 
of back-reaction is precisely what we expect by looking at equation (p7| ), due to the presence of the 6G^ term. When 
the energy density perturlaations of the total fluid are dominated By perturbations in the matter component (i.e . 
at sufficiently high redshifts in typical dark energy cosmologies) for some scale fc~^, the term 5G^ — k6Tq in ( |34| ) 
is in turn dominated by matter energy density perturbations, which then act as a source of the scalar field density 
perturbations. 

The very interesting feature here is that non-vanishing energy density perturbations of a non-minimally coupled 
scalar field can even be associated to a homogeneous scalar field as long as a non-zero value of (j) makes k^^ 7^ F: we 
can easily see that Spff, in equation (^) survives even in the limit S4> —> 0, because of the gravitational dragging. In 
other words, perturbations of a non-minimally coupled scalar field are sourced by two complementary mechanisms: 
proper scalar field perturbations, and metric induced perturbations, related to the Ricci scalar coupled with the field 
itself. 

Focusing now on dark energy cosmologies, the described process introduces genuine new features with respect to 
"ordinary" (minimally coupled) Quintessence scenarios: the growth in the matter perturbations may drag EQ density 
perturbations to a non-linear regime, opening the possibility of the formation of Quintessence clurnps. 

To give a concrete example, we numerically evolve linear perturbations in a typical EQ scenarios M. The coupling 
of Quintessence with the Ricci scalar is chosen to have the structure 

-fi^,R) = Fi^)R. (41) 

K 

The measured gravitational constant G. in scalar-tensor theories (0) with the choice (^l]), is related to the various 
quantities in the Lagrangian as follows [31 : 



kF V 2ujF + 3F^^ 

As in ||r^ , |20|] , we model _F as a constant plus a term yielding a quadratic coupling between the field and the Ricci 
scalar, so tmit 

F(0) = -+e0^ (43) 

K 



where ^ is the non-minimal coupling constant, with the constraint that today F satisfy the relation (|42|). Note 
also that ^ can in principle assume both positive and negative values. Moreover, its magnitude is not arbitrary, 
due to constraints from local and solar-system experiments on the time-variation of the gravitational constant and 
from effects induced on photon trajectories j3^]. The time derivative Gt of the observed gravitational constant G as 
defined in equation ( p2[ ) must satisfy Gt/G ^ 10^^^ yr^^ at present. Moreover the Jordan-Brans-Dicke parameter 
wjBU — ^F/F'^^ must be greater than about 2500 at present; to be conservative, we choose lojbd = 3000, with a 

negative sign of the coupling constant (in our specific model this corresponds to ^ ~ —1.78 • 10~^, 0o — I/a/G), which 
yields Gt/G ~ 10"^** yr~^. The Quintessence potential, responsible for cosmic acceleration today, has an inverse 
power law form V cx moreover, we fix Lij{4i) = 1. 

The cosmological model is specified as follows. The Hubble parameter at present is fixed at Hq — lOO/i km/sec/Mpc 
with h = 0.7 and the spatial metric is taken to be flat, K = Q. The fraction of critical density in dark energy is 
= 0.70. The equation of state of Quintessence at present w^o is chosen to be —0.9, yielding cosmic acceleration. 
Baryon abundance is set to Q.bh^ — 0.022, CDM represents the remaining matter component, ^cdm = ^ — and 
three massless neutrino families are assumed. Perturbations are taken to be Gaussian with an initially scale invariant 
adiabatic spectrum ||37|] . The evolution of background and perturbations has been determined by numerically solving 
equations in the synchronous gauge A = B ~ 0. Their expressions are reported in the Appendix. 



8 



Let us consider the background evolution first. In figure I the energy densities of radiation (dashed hue), matter 
(dotted) and dark energy (sohd) are plotted as a function of the redshift z = 1/a — 1. At late times, z ^ 5, the dark 
energy density is dominated by the kinetic and potential energies. At higher redshift the effect of the gravitational 
dragging is evident: the last term in (O) actually dominates and forces dark energy to scale with redshift as the 
dominant cosmological component. As it is easy to see, in this regime, the dark energy cosmological abundance is 
simply given by 

n^z) ~ -kC^^z) , (44) 

where the minus is due to our sign conventions. Note that the quantity can be constrained by Big Bang Nucle- 
osynthesis (BBN), because a variation in the gravitational constant can be regarded as inducing a change into the 
effective number of massless neutrinos [p8[ ; however in the present case this value is at percent level, too small to 
produce observable effects. The dotted-dashed line represents indeed only the contribution from the mc terms in ([l^): 
the rising part of this curve at z ^ 1000 is due to the _R-boost induced by the effective gravitational potential in 
the Klein-Gordon equation (|l8|). We stress that while the latter contribution comes from the kinetic energy of the 
field rolling on the effective gravitational potential, so that ultimately implies a change in time of the physical value of 
the scalar field </>, the gravitational dragging can be thought as a power injection into the dark energy density coming 
from the total one, while it does not require directly a spatial dependence of the expectation value of (/). Note also that 
a condition in which dark energy scales as the dominant cosmological component can be achieved with an exponential 
potential ]^,[Tot; however, in that case the field is minimally coupled and Quintessence density perturbations vanish 
relativistically after horizon crossing p^ ]. 

Let us turn now to consider the perturbations. Since we are interested in the dark energy clustering during the 
formation of matter structures, we concentrate on the behavior in the matter dominated era. Moreover, we focus on 
the logarithmic power of density fluctuations at the scale k, defined by 

Sl^^Ank^Sl, (45) 

where x represents a generic component. As it is well known, a scale for which Skx — 1 has to be considered in 
non-linear regime. In figure ^ we plot 6k for matter (dotted) and dark energy (solid) at relevant redshifts. As it 
is expected, the gravitational dragging is active and forces Quintessence perturbations to behave as non-relativistic 
matter on sub-horizon scales, when the gravitational terms dominate both and Sp^. Under these conditions, by 
using ( |l6|) and (|3^) we see that we can write approximately 

S4,^S^, (46) 

which is mostly satisfied, in our specific model, at 2 ^ 5. It is in fact not a case that all the quantities specifying 
our specific model disappeared in the above relation. Indeed, (^) holds if three general conditions are satisfied in 
dark energy cosmologies, namely (i) Gravity deviates from General Relativity, (m) Quintessence plays the role of 
the non-minimally coupled field and (Hi) gravitational terms dominate (^,|l6|). At present, when the kinetic and 
potential energies of the field dominate the cosmic expansion imprinting acceleration, the condition ( ^ ) is broken 
because matter is no more the dominant component. 

To fully understand the importance of this plot, we reported also Skip for a corresponding model in which the 
Quintessence is minimally coupled, represented by the dashed dotted curve. The latter is rising with time because, as 
noticed in p9| ] , the inhomogeneous term of the perturbed Klein Gordon equation ( p^ ) is driving the evolution of S(f>. 
However, as it is evident from the figure, in absence of non-minimal coupling the dark energy density perturbations 
do not play any role in structure formation. 

Therefore, maybe the most interesting consequence of the gravitational dragging in dark energy cosmologies is that 
the non-linearity may arise for the Quintessence component, at a redshift depending in particular on the coupling 
strength, opening the possibility of the formation of Quintessence large overdensities and cavities on sub-horizon 
scales. On the other hand, at a linear level, the effect produced by on the total gravitational potential <i> resides in 
a fraction p,p/pm, which is small in the limit in which is subdominant. For example, in models in which the assumed 
scalar-tensor Gravity theory is slightly different from General Relativity, i.e. |1 — K^Fj 1, in the gravitational dragging 
regime when the gravitational contribution dominate and through the products {k~^ — F)Gq and (k^^ — F)SGq, 
it can be easily seen that the the portion of gravitational potential which is sourced by the Quintessence is given 
approximately by 

$0 ~ (1 - Ki^) $ < $ , (47) 

which means that the bulk of the gravitational potential arising from the clustering process is still provided by matter. 
Note however that this is true for linear perturbations, which we are treating in this work; the effect a non-minimally 
coupled dark energy on the gravitational potential associated to a non-linear structure is still unknown, and actually 
this study could be an interesting development of the present work. 
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The possibility of the presence of Quintessence clumps in the Universe has been considered by several authors, 
mainly aimed to foresee their signatures on the Galactic structure p3|-p5[|; however, there is not, at the present, a 
theory explaining how such "vacuum energy clumps" may form, and on wmch scales we expect eventually to find them 
at low redshifts. What we obtained is just a possible way to escape from the linear regime. In other models, as in 
p7| , there can be even higher deviations from General Relativity at high redshift, than in the model considered here; 
fHe potential effect on the formation of non- linear clumps may be even more important. Moreover, we have shown 
that the occurrence of dark energy non-linear sub-horizon structure does not require directly a space dependence of 
the expectation value of (j). 

Obviously, this is just a first step towards a theory of "vacuum energy clumps" : the perturbation behavior out 
of the linear regime is still an open issue, and is strictly related to properties such as the sound speed of the dark 
energy component. In particular, we think that one of the key issues is the time needed for a Quintessence primordial 
structure to collapse, and a fundamental role on this is played not only by the bare coupling constant ^, but generally 
by the whole coupling function, through the effect they have on the value, and sign, of the Quintessence sound speed. 
In the next S ection we'll use a different approach to explain such properties, based on the effective sound speed 
introduced in |2§| , |34|| . 



V. DARK ENERGY SOUND SPEED 



An important role in the perturbation growth is played by the non adiabatic stress or entropy contribution, entering 
directly into the evolution equation for density perturbations; by defining Spx — PxT^Lx for a generic component x we 
can write the scalar field entropy contribution as 

p^T^ = 6p^ - clJp^ , (48) 



where i is the adiabatic sound speed of the scalar field, defined in equation (22), which can also be written as 



4 = ^^~i7T^- (49) 

In most Quintessence scenarios, the field is modeled as a component with negative, and slowly varying, equation of 
state, so that from eq. (^, ~ w^. As discussed in p8| , ^ , looking at the continuity and Euler equations ( p^pO[ ), 
a fluid with negative sound speed without entropy and stress terms would make adiabatic fluctuations unable to give 
a pressure support against gravitational collapse of density perturbations: in other words, under these circumstances, 
the adiabatic pressure fluctuations would accelerate the collapse rather than oppose it, as can be derived from the 
density perturbations evolution in the limit of sub-horizon scales. In this scenario, density perturbations would rapidly 
become non-linear after entering the horizon, unless the cntropic term in eq.(p^) acts as a stabilizing mechanism: this 
requires w^F^ > 0. To check this possibility, W. Hu in [|8| introduces the effective sound speed c^yj ^ defined in the 
rest frame of the scalar field, where ST^^ = (in |2^, these considerations are applied to a more general "generalized 

dark matter" component, which can recover the Quintessence scalar field case, as well as matter and radiation). 
Following this approach, we write the gauge invariant entropic term as 

z.,F, = (c2,,-c2)5('--*) , (50) 

where J^'^'^**-' is the density contrast in the dark energy rest frame, which therefore coincides with the gauge invariant 
density perturbation as follows: 

Sl;'"'^ =A^ = S^ + 3jil + w^) [v^ - B) . (51) 

By doing so, the stabilization mechanism of scalar field perturbations is expressed only in terms of gau ge invariant 
quantities. In the mentioned case of c^^ < 0, effective pressure support is obtained if the entropic term (|50| ) is positive. 
The effective sound speed can be interpreted as a rest frame sound speed; importantly, it allows to define a stabilization 
scale for a perturbation, given by the corresponding effective sound horizon. This formalism has been used in pq ] 
to show that density perturbations in a minimally coupled scalar field of Quintessence are damped out below tne 
horizon, so that the Quintessence rapidly becomes a smooth component: in this case, indeed, it can be verified that 
the effective sound speed is c^ff^ = Ij giving a relativistic behavior to the corresponding density fluctuations. 

The situation in Extended Quintessence scenarios can however be much different, in what the effective sound speed 
may be strongly affected by the presence of the non-minimal coupling term. Rewriting c^^^ as 
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2 _ 6p^ + cl^mh^{v4,~ B)lk 



where /i^ = Pcji, it is quite evident that, on sub-horizon scales fc 3> 7i, c^^j^ ~ Spcj,/Sp^. As discussed in the 
previous Section, in dp^ lies the main difference between ordinary Quintessence and Extended Quintessence: while 
5pij)/5p^ ~ 1 for minimally coupled scalar fields, giving relativistic values to clj^j and damping field perturbations, 

this ratio may be much lower than unity whenever the energy density perturbations of the scalar field are enhanced 
by perturbations in the matter field, and this is a peculiar property of non-minimally coupled scalar fields. 

Another genuine feature which is expected in EQ scenarios regards the viscosity of the dark energy component. 
As pointed out, again in |p8[|, the anisotropic stress can also be a smoothing mechanism for the scalar field, damping 
density perturbations through its effects on velocity perturbations in equation (po|). A viscosity parameter c^^^ is 
introduced to relate velocity /metric shear and anisotropic stress; for the Quintessence scalar field, we have 



4c^ 

ttt + STYttt = — "-^[kvd,-HT] ■ (53) 
W0 V / 



In the limit of negligible ttt cj, and in shear- free frames {Ht — 0) , c^j^ > determines a viscous damping of velocity 
perturbations, as it can be seen through the Euler equation (^), which sums up with the viscosity effect arising from 
the cosmological expansion; thus, if c^^^ > 0, viscosity can be an extra smoothing mechanism. On the other hand, 
c^is < results in a term which acts against the cosmological viscosity into the Euler equation. 

The viscosity parameter turns out to be zero for minimally coupled Quintessence, where anisotropic stress is not 
present: in that case, however, c^jy = 1, so that the adiabatic stress only is enough in smoothing the scalar field on 
sub- horizon scales. On the other hand, for Extended Quintessence fields, we expect a non zero contribution to c^^^, 
due to the traceless part of (5T;"™=[0] -I- JT]»™"[0]. 

We plot the four quantities Cg^-j^, w^, c^^, c^j^^ in figure^, comparing results in our tracking EQ scenario (solid 
lines) and in an ordinary minimally coupled Quintessence cosmology (dotted dashed). 

The most striking effect is for c^j^^. For all the redshifts relevant for structure formation the effective dark energy 
sound speed is vanishing in the EQ case, allowing for a behavior of its density perturbations analogous to that of 
non-relativistic matter. Correspondingly, minimally coupled Quintessence has c\ff^ — 1- This reproduces the same 
result as in figure ^ obtained with a different approach. The more the gravitational term in 5p^ dominates in the 
denominator of the expression (p3), the larger is the suppression of the dark energy pressure reaction to the density 
contrast growth. Even if the plocted results are strictly valid only in our model, we stress that this is an example 
of a general occurrence in scalar-tensor dark energy cosmologies. In addition, since the gravitational term in ( p7| ) is 
proportional to {\/ k — F) which can in general assume both positive and negative values, even the sign of c^ff^ can 
be reversed realizing a scenario in which the sound speed accelerates the collapse on scales smaller than the sound 
horizon instead of opposing it. 

An analogous behavior can be found by looking at the plots of and c^^, in what that they are severely depressed 
at relevant redshifts in EQ models. Also, the difference between c^^ and is due to the time derivative of the equation 
of state through equation (|49|). Notice in particular that in minimally coupled Quintessence we have cl^,w^ < 
while Ceffrf, = 1; as we already stressed, the latter quantity is indeed the appropriate one to explain the behavior 
of a minimally coupled scalar field, resulting in the relativistic damping of sub-horizon perturbations. Finally it is 
interesting to check whether viscosity can be sufficiently effective in damping out density and velocity perturbations, 
even when adiabatic pressure fluctuations are not. As in the minimally-coupled case, it turns out that this is not the 
case; first, the amplitude of c^^^ ^ is much lower than unity, and second the negative sign yield an enhancement of 
velocity perturbations, instead of a damping, as we stressed above; the peak at z ~ 1 is due to the onset of cosmic 
acceleration, i.e. the sign change into the cosmic equation of state. 

In practice, there are no mechanisms to slow or decrease the amplitude of Extended Quintessence density perturba- 
tions in the gravitational dragging regime. This analysis confirms, and better clarifies, the results of previous Section 
concerning the scalar field density fluctuations power spectrum. 
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FIG. 2. Spectral power of density fluctuations for matter (dotted) and dark energy in tracking EQ (solid) and Q (dotted 
dashed) scenarios at different redshifts. 



VI. CONCLUDING REMARKS 



We studied the behavior of linear perturbations in scalar-tensor cosmologies, focusing on the density fluctuations 
of the scalar field (f> coupled with the Ricci scalar R. We found that such coupling can activate a "gravitational 
dragging" of the scalar field density fluctuations by the cosmological metric perturbations, which in turn are powered 
by the whole cosmological stress-energy tensor through the Einstein equations. That is, as the non-minimal coupling 
represents a power exchange channel between the scalar field component and the general relativistic cosmological 
gravitational potentials, we studied in particular how such channel acts at the level of linear density perturbations 
in the scalar field, represented in particular by the density contrast = 5p^l pcf,. In conditions in which (p is not 
the dominant cosmological component, the power injection coming from Gravity can largely dominate 5^ forcing its 
dynamics to be similar to that of the dominant component itself. On the other hand, in the same conditions, the 
scalar field contributes to the cosmological gravitational potentials by a fraction given by the ratio between the scalar 
field and total energy densities. 

This phenomenology has important consequences on the dark energy clustering properties in Extended Quintessence 
scenarios, where the non-minimally coupled scalar field is assumed to be responsible for the cosmic acceleration today. 
Namely, the dark energy assumes the features of a pressureless fluid when non-relativistic matter m dominates, i.e. 
after matter radiation equality and in the pre-accelerating stage of the cosmic expansion. In other words, the scalar 
field density perturbations can grow on sub-horizon scales, tracing those in the matter component; this fact is depicted 
in the scalar field density contrast (5^, as well as in the properties of its effective sound speed Ce//^: 

5^~5.,n , Ce//0<1 . (54) 
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As we already mentioned, the reason of this behavior Ues in the gravitational coupling to the Ricci scalar, contributing 
a gravitational term in the scalar field energy density which gets the dominant contribution from the perturbation in 
the matter component. The latter perturbations are therefore able to feed the dark energy density fluctuations up to a 
large amount even if the non-minimal coupling is small enough to respect all the existing constraints on scalar-tensor 
theories of Gravity. We stress also that the behavior (p4) is not depending on the particular form assumed to describe 
the non-minimal coupling; indeed, such gravitationardragging regime holds whenever the contributions due to the 
non-minimal coupling dominate both in and Sp^, so that their ratio is rather insensitive to the detailed shape 
of such coupling. Moreover, it should be noticed that the behavior (blh is not to be related to variations of the 
expectation value of the scalar field 0; indeed, our study shows that densTEy perturbations of a non-minimally coupled 
scalar field are sourced both by fluctuations of expectation value as well and by perturbations of the Ricci tensor. In 
particular, the dragging effect emphasized here is generated even in the limiting case of a homogeneous scalar fleld, 
being induced by the coupling with R. 

We have provided a worked example of the above phenomenology in Extended Quintessence scenario, involving a 
quadratic coupling between the field and R. Numerical integrations of the cosmological equation system shows that 
the dynamical condition is satisfied at redshifts relevant for the structure formation process, respecting all the 
existing constraints on scalar-tensor Gravity theories. 

We believe that these results open new perspectives on the standard picture of structure formation in dark energy 
cosmologies, since the gravitational dragging expressed by (j^) implies that both dark energy and matter exit the 
linear regime on sub-horizon cosmological scales at relevant redshifts. This immediately poses the problem of their 
evolution afterwards, i.e. the gravitational clustering of large overdensities and deep cavities composed by matter 
and scalar energy tangled by a non-minimal gravitational interaction; while as we already stressed the gravitational 
dragging is rather insensitive to the detailed shape of the non-minimal coupling, the same could be untrue at a non- 
linear level. In particular, for a given model, it would be interesting to look at the appearance of the resulting density 
profile after virialization, since this aspect could be constrained by observed rotational curves in nearby galaxies. 
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FIG. 3. Redshift behavior of effective sound speed, equation of state, sound speed and viscosity in tracking EQ (solid) and 
Q (dotted dashed) scenarios. 



VII. APPENDIX: SCALAR FIELD PERTURBATIONS IN SYNCHRONOUS GAUGE 



Numerical integration and analytic treatment of the perturbation equations get simplified when developed in the 
synchronous gauge pO{ . To obtain our numerical results, we used a modified version of cmbfast which covers 
Extended Quintessence cosmologies. 

A scalar-type metric perturbation in the synchronous gauge is parameterized as 



(55) 



hij{x,T) 



kikjft,(k, r) + (kjkj - -(5ij)67y(k, r) 



(56) 



being h the trace of hij . By choosing A = B = in ( [T2| ) , the metric perturbations and Ht are related to h and 
77 by the following relations: 



2 ' ' 6 



(57) 



Defining the shear perturbation cr by the relation a = ai/k^piTTY respectively, one has, for the scalar field, the 
following quantities derived from the conserved T^flip] (eq. y), in synchronous gauge: 



6p = 5p fluid + + o 



1 / 



K 



nSF 



R 



Fh 
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6a'' \ K 



-nh-- + k^T] 



(59) 



{p + p)v = {p fluid + Pfluid)v fluid + 



oj(t)64> + 5F- H6F + 2 



F]r, 



(60) 



where 



(P + P)<^ = {P fluid + Pfluid)o-fluid + 



+ 



3a2 



1 rA I '^i 6W . 2 •■ 12 .. 



P^Ti'L is the isotropic pressure perturbation. The perturbed Klein-Gordon equation reads 

SR . 



a; / ,0 2 



2cj 



2cj A: 



These perturbations enter in the perturbed Einstein equations, easy to solve in this gauge: 



k^r]- -nil = - 



/.2^ _ a'<P + P)v 



h + 2nh - 2k^T] = -3a^K6p , 
/i. + 6f/ + 2n{h + 6f)) - 2k'^r] = -Sa^Kip + p)a 



(61) 

(62) 

(63) 

(64) 

(65) 
(66) 



This set of differential equations requires initial conditions on the metric and fluid perturbations; we adopt adiabatic 
initial conditions ||37[|). 
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